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We introduce and study systems of randomly coupled maps (RCM) where the relevant parameter is the degree of connectivity 
in the system. Global (almost-) synchronized states are found (equivalent to the synchronization observed in globally coupled 
maps) until a certain critical threshold for the connectivity is reached. We further show that not only the average connectivity, 
but also the architecture of the couplings is responsible for the cluster structure observed. We analyse the different phases of 
the system and use various correlation measures in order to detect ordered non-synchronized states. Finally, it is shown that 
the system displays a dynamical hierarchical clustering which allows the definition of emerging graphs. 



PACS number(s): 05.45.+b, 05.20.-y, 05.90.+m 
I. INTRODUCTION 

Since their introduction in 1989 fj],||], globally coupled 
maps (GCM) have turned out to be a paradigmatic exam- 
ple in the study of the emergent behavior of complex sys- 
tems as diverse as ecological networks, the immune sys- 
tem, or neural and cellular networks. It has been shown 
that mutual synchronization of chaotic oscillations is a 
robust property displayed by GCM (|||. Effects of mu- 
tual synchronization are also known for coupled chaotic 
oscillators with continuous time 0. A recent study || 
of large globally coupled populations formed by chaotic 
Rossler oscillators has revealed that mutual synchroniza- 
tion and dynamical clustering in these systems are simi- 
lar to the respective behaviour found in GCM. Another 
well investigated class of self-synchronizing systems rep- 
resents lattices of locally (e.g., diffusively) coupled oscil- 
lators Moreover, oscillator systems with both local 
and global coupling between elements have also been dis- 
cussed . A common property of GCM and other above- 
mentioned systems is their architectural symmetry: The 
pattern of connections of any individual element is iden- 
tical. This internal symmetry is preserved in the fully 
synchronized dynamical states and spontaneously broken 
when dynamical clustering takes place. 

The architecture of dynamical networks found in real 
ecological or other systems would rarely be so sym- 
metric. Even in the situations with high connectivity, 
when the links can extend to many distant elements, the 
graph of connections may have a complex topology. The 
question is therefore whether and in what form mutual 
synchronization and dynamical clustering can persist in 
such complexely connected networks, lacking a structural 
symmetry. In this direction, among the few cases already 
explored we can mention a layered system of identical 
random neural networks with partial (though regular) 
connectivity among layers || , coupled lattice maps with 
connections extending further than to nearest neighbors 
M, as well as coupled lattice oscillators in 2-D with dif- 
ferent coupling schemes ]Tq] or a model ecosystem with 



partial connectivity ETJ. More recently, ensembles of 
nonlinear oscillators [ 12]] and GCM |13| with random in- 
teractions and variable symmetry (still globally coupled) 
have also been analyzed. 

The analysis of systems formed by few coupled caotic 
elements complements the previous approach and has 
provided a better understanding of the collective be- 
haviour displayed by large ensembles of regularly con- 
nected chaotic oscillators. The stability of the syn- 
chronous state can be in this case easily quantified by 
means of the transverse Lyapunov exponent Aj_, which 
changes sign in a blow-out bifurcacion |14j] and makes 
the synchronous state unstable. Near the bifurcation, 
on-off intermittency fTq ] and riddled basins of attraction 
are observed. Much attention has been devoted to 
the dynamical properties of systems formed by two cou- 
pled, identical chaotic elements \LT^ . The synchroniza- 
tion properties of two coupled logistic maps have been 
extensively investigated JT8| . 

As the next step towards understanding the synchro- 
nization phenomena in complex networks, one can con- 
sider systems formed by a large number of identical dy- 
namical elements that are connected by identical sym- 
metrical links but where the pattern of connections be- 
tween elements is random. This is the starting point 
of our paper. In the following sections we introduce 
and analyse what we have termed Randomly Coupled 
Maps (RCM), that is, networks of chaotic maps con- 
nected (symmetrically) at random where the relevant pa- 
rameter is the average connectivity in the system. 

Our main result is that mutual synchronization and 
dynamical clustering are possible in RCM even when a 
significant fraction (up to 40-45%) of all potential connec- 
tions is absent. However, the synchronization and clus- 
tering phenomena in these systems are different in certain 
aspects from what is known for GCM. Exact synchro- 
nization and the formation of identical dynamical states 
of the elements are not found here. Instead, either one 
or several compact clouds (fuzzy clusters) are formed. 
These fuzzy clusters are dynamical objects which split 
into subclusters or join other groups of elements. Such 
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dynamical hierarchy of clusters is almost never com- 
pletely fixed in time. A closely related effect is that the 
asymptotic dynamical behavior in RCM is never sensitive 
to the initial conditions. However, the architecture of a 
particular network may bias the synchronization process 
and make certain cluster distributions more favourable. 
The role of the network is particularly evident when the 
system is small, (N ~ 0(1)): In this case, the synchro- 
nization properties of the system are strongly dependent 
on the particular architecture, and graphs with the same 
connectivity might have very different collective behav- 
ior. We will show that only in the thermodynamic limit 
do the synchronization properties of RCM become equiv- 
alent to the globally coupled case. 

In the next Section we introduce the model and de- 
scribe its dynamical behavior. In Section III the syn- 
chronous and the partially ordered phases of the system 
are quantitatively characterized by computing properties 
such as the mutual information and two order parameters 
of the synchronization transition. A more detailed statis- 
tical investigation is then performed in Section IV where 
distributions over pair distances are constructed. Sec- 
tion V is devoted to the analysis of the emergent cluster 
structure in partially condensed phases. A complement 
to the latter section is the Appendix, where an example 
of a small system with varying network architectures is 
considered. Finally, in Section VI we discuss our results, 
come to the conclusions and outline future extensions of 
this study. 



II. THE RANDOMLY COUPLED MAPS 

We begin by exactly defining what we have termed 
Randomly Coupled Maps. Instead of a globally coupled 
system, this will be a network of connections character- 
ized through a random matrix Ty , the elements of which 
are either (when a connection between maps i and j is 
absent) or 1. In our analysis, we assume that the matrix 
is symmetric, i.e. Ty = Tji, and all diagonal elements 
are set to (Tn = 0). An important property of such 
random networks is their average connectivity 
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Hence, each element will be on the average connected to 
v(N — 1) elements in the system. If v = 1, the system is 
globally coupled and our system reduces to this known 
system. 

The collective dynamics of RCM is defined as 
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where e specifies the strength of the coupling and f{x) 
is the individual map. This collective dynamics can thus 
be understood as involving diffusion on a graph: Each 
node in the random graph diffuses a certain fraction of 
material (its state) to the elements to which it is con- 
nected, and receives a comparable influence from them. 
The intensity of coupling between the network elements 
is specified by the parameter e (we have the constraint 
e < v). Nevertheless, the combination e/v, which can 
be interpreted as a diffusion coefficient, plays an impor- 
tant role. Note that in the limit vN — ► oo all RCM 
with a given connectivity v become statistically identi- 
cal. Indeed, in this limit each element will have the same 
number vN of connections, although still it is linked only 
to a randomly chosen subpopulation. If the mean-field 
approximation holds, the behavior of this system would 
therefore be equivalent to the behavior of a GCM with 
coupling intensity e. 

In this paper we investigate differences in the behavior 
of GCM and RCM with large, but finite numbers of ele- 
ments. We work with the logistic map f(x) = 1 — ax 2 and 
use values of a such that the dynamics of a single map 
is chaotic. Our numerical simulations are performed for 
graphs with sizes up to N = 2048 elements and with con- 
nectivities in the interval 0.5 < v < 1. Such graphs have 
been randomly generated by independently choosing ev- 
ery possible connection with the same fixed probability. 
We have checked that thus generated graphs remained 
fully connected, i.e. they could not be further separated 
into two disconnected parts. 

Considering that we wish to compare the behaviour 
of RCM with that of GCM, it is of interest to begin by 
briefly recalling the dynamic behaviour of GCM when 
the coupling strength e is decreased from 1 to (cf. (p|). 
When, for example, a = 2, the synchronous phase in 
GCM is maintained until e = 0.5, where it destabilizes 
and is substituted by the so called 'glassy phase'. In 
Kaneko's terminology, 'glassy' means that the final at- 
tracted is sensitive to the initial state of the maps, and 
thus a multiplicity of attractors is to be found in this 
phase for the same parameters value. For e w 0.33, an 
ordered phase sets in: The sensitivity to the initial condi- 
tions disappears and the elements group together in only 
2 or 3 different clusters. A narrow intermittent band in 
0.21 < e < 0.25 preceeds the turbulent phase, in which 
the number of clusters in all the attractors is of order 
0(iV). 

We first examine the possibility of full mutual synchro- 
nization in RCM. Such full synchronization takes place 
in GCM when the coupling strength exceeds a critical 
value ||. In the fully synchronized (coherent) regime the 
states of all maps in GCM are identical. Our analysis 
reveals that such exact synchronization does not occur 
in RCM. However, at sufficiently high coupling intensi- 
ties, all elements move together in a single compact cloud 
(a fuzzy cluster), so that the typical distances between 
their trajectories remain below \x l (t) — x J (t)\ < 10~ 8 
(this is the single precision for real numbers in our com- 
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puter). The fuzzy synchronous phase is maintained for 
a certain range of values of the coupling strength e, and 
then a sudden transition to an asynchronous state (see 
the discussion below) is observed. Fig. 1 shows how 
this transition proceeds for a randomly chosen network 
with v — 0.8 and N = 50 elements. For comparison, 
we show in the same figure the respective behaviour 
of GCM. While for GCM the turbulent-ordered-glassy- 
synchronous sequence of phases is clearly seen when e 
decreases, only the synchronization breakdown for RCM 
is apparent. 




FIG. 1. Numbers of clusters as functions of the coupling 
strength e for RCM of size N — 50 and connectivity v — 0.8 
and for GCM (y = 1) of size N = 50. In the RCM case, two 
elements i and j are considered to belong to the same cluster 
if Ix'ft) - x ] {t)\ < 10~ 8 . In the GCM case, two elements 
belong to the same cluster if their states are identical. 

We have studied how the synchronization threshold in 
RCM depends on the mean connectivity v of the network 
and on its size N. We accept that the elements in a 
given network are synchronized if the average distance d 
between pairs of elements during an interval At is smaller 
than 10 -8 , after a transient is discarded, that is 
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The simulations begin at a value of e < e* which is step- 
wise increased in amounts 10 -3 until the condition (|^) is 
fulfilled. We have also observed that immediately prior 
to the synchronization transition, the distance d under- 
goes a sharp cutoff and afterwards stabilizes around small 
values which depend on the system size (d is in the syn- 
chronous phase a decreasing function of N). In Fig. 2 
the distance to the synchronization threshold of GCM is 
shown as a function of the network size when the con- 
nectivity is kept constant (y — 0.9). Numerically, we 
obtain 



Large open circles are averages over 10 to 30 independent 
graphs, and are represented together with the dispersion 
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in the thresholds (error bars). The index j stands for dif- 
ferent networks formed by the same number of elements 
N, while e*(N) is the average value for each size. Small 
solid circles represent the synchronization threshold for 
each of these graphs. The error in the determination of 
e* for each network is 10~ 3 . The solid line in the main 
plot has slope —1/2. 
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FIG. 2. Dependence of the synchronization threshold e* on 
the network size N. The open circles correspond to aver- 
age values of the critical coupling intensity e* for networks of 
sizes ranging from N — 2 3 to 2 11 with the same connectivity 
v — 0.9. The solid line has slope —1/2. In the insert we show 
the dispersion in the values of e* as a function of the system 
size. The solid line has slope —1. In this case, a — 2 and thus 

e *GCM = 0.5. 

As larger networks are considered, the dispersion 
Ae*(N) in the synchronization thresholds for different 
networks with the same connectivity becomes smaller 
(insert). Our numerical results point to a dependence 
of the form Ae* ~ iV _1 . The fact that the dispersion in 
the values of e* tends to zero in the limit N — > oo for 
a fixed v indicates that RCM are characterized by self- 
averaging quantities. In view of these results, the RCM 
should be well described by a mean-field approximation, 
thus by GCM, in the limit vN — > oo. In fact, the numer- 
ical simulations agree with this picture. In the opposite 
limit, when the number of elements in the systems is 
small (N ~ 0(1)), the threshold at which the group first 
synchronizes is very sensitive to the particular way the 
maps are connected, as can be already seen for TV < 64 
in Fig. 2 (see also the Appendix). Fig 3. depicts the 
synchronization threshold for systems with N — 50 and 
N = 500. Every point corresponds to a fixed network 
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with connectivity given in the x— axes and the critical 
value e* in the y— axes. 
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FIG. 3. Threshold to synchronization. Small solid circles 
are critical thresholds for systems of size N = 50 and big open 
ones correspond to N = 500. We restrict ourselves to values 
e < v, as discussed in the main text. The line e = 1/2 is the 
stability threshold for GCM, and acts as a lower boundary for 
RCM. 
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FIG. 4. Divergence in the transient time t c when approach- 
ing the critical point e* for a randomly chosen network with 
TV = 50 and v — 0.8 (squares) and for the GCM of the same 
size (filled circles). The solid line indicates the divergence 
law with t — 1. The inset displays the number of clusters as 
function of time for the RCM at e = 0.7. 

The synchronization threshold for GCM has previously 
been exactly determined ||. When N — > oo it is given 
by the condition e = £qq M where the critical coupling 
is obtained from A + log(l — ^gcm) < ano ^ ^ ^ s the 
Lyapunov exponent of the single logistic map. It can 
be expected that the synchronization threshold for RCM 
approaches the limit e = £* G cm wnen v — ► 1, for any 
value of N, and also when vN — ► oo, and independently 
of v : as has been discussed. It can be seen from the 
numerical results represented in Figs. 2 and 3 that the 



value e = 0.5 (for a = 2) yields indeed a lower estimate 
for the synchronization threshold in RCM. 

The system falls into the synchronous phase after a 
transient of diverging length when approaching the syn- 
chronization threshold e*. Fig. 4 shows the dependence 
of the transient time t c on the distance e — e* to that 
threshold. The insert shows the dynamics during a typi- 
cal transient. We see that strong fluctuations (intermit- 
tency) are accompanying the convergence process. Rep- 
resenting the dependence of the transient time in the 
form t c oc (e — e*)~ T , we find that the exponent r sa 1 is 
typical both for RCM and GCM. More precisely, a least 
squares fit to numerical data returns tqcm = 0.997(3) 
and trcm — 0.982(3), and the critical threshold values 
to synchronization are £qcm = 0.5 and e* RCM /u w 0.87 
for the particular graph of Fig. 4. 

Below the synchronization transition, the glassy phase 
is observed in GCM. The dynamical behavior of RCM in 
the region e < e* is investigated in the next section. 

III. MUTUAL INFORMATION AND THE TWO 
ORDER PARAMETERS 

We have used three different measures of correlations 
among elements to check if the phase with e < e* has 
still some intrinsic order . The first of them is the mutual 
information between two maps i and j, P 3 (£). To define 
this quantity, we introduce a partition of the phase space 
of the logistic map in the following way. If the state of 
the chosen element i is x % (t) > then it will be assigned 
a value 1, and if x l (t) < 0. This generates a sequence 
of bits in a certain time interval SI € £ = {0, 1} which 
allows the calculation of the Boltzmann entropy for the 
ith map ]l9] , 

= - ]T P(Sl)\nP(Sl). (5) 

S|=0,1 

In a similar way we define the joint entropy for each pair 
of maps, 

= - J2 E P(Si,Si,)lnP(St,Si,) (6) 

S t l =0,l S j . =0,1 

and finally the mutual information for i and j is given by 
I ij (S) = + H j (Y.) - H ij (Y,). (7) 

The mutual information is a good measure of correla- 
tions, e.g. it achieves maximal values near critical points 
Ipof . In a context closer to ours it has been shown to ac- 
curately discriminate among the different phases of GCM 
f2lf . An interesting property of P 1 is that it is practically 
precision-independent, due to the rough coarse-graining 
of the dynamics. 

As an illustration of the sensitivity of this measure, 
let us discuss which values of P 3 are expected in two 
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dynamically opposite regimes, i.e for the synchronous 
and the turbulent phases. In the coherent (synchronous) 
phase where the states of all elements are identical, the 
sequences SI are the same for all of them, therefore 
if'(£) = £P(£). In addition, since the two chosen maps 
are visiting the same points, we have = H l and 
thus Pi = H l , reflecting the trivial nature of the cor- 
relations. When a — 2, because of the symmetry in the 
invariant measure of the logistic map for this parame- 
ter (and noting that the synchronized system is equiva- 
lent to the single map) the mutual information achieves 
its maximum value. In fact, under these conditions, 
F(0) = P(l) = 1/2, thus maximizing W = ln2 and 
also P 3 ' — In 2. For parameter values a ^ 2 the invari- 
ant measure is not symmetric around 0, -P(O) ^ -P(l) 
and hence P 3 < In 2 typically. Nonetheless, as can be 
seen in Fig. 5, the synchronous state is clearly detected 
through this measure. In the turbulent phase (assuming 
that the elements behave independently and are not cor- 
related) we will have again H 1 (Yj) = £F(£), but now the 
joint probabilities factorize, P(Si,S 3 ,) = P{Sl)P(S 3 ,) 
and TP 3 = 2H l . In this phase we therefore expect P 3 « 
irrespectively of the parameter a. In the intermediate 
cases, where some correlations are present, the mutual 
information should take values between the former two 
limits, < P 3 < In 2. 

To determine the mean mutual information 

(pi) = —Jl V" pi ( 8 ) 

1 ; <*<i> 

for a network with a given matrix Tij , we take an average 
of the mutual informations for all possible pairs (i,j) for 
long enough sequences to ensure the stability of the prob- 
abilities (typically At = 10 3 — 10 4 after discarding a tran- 
sient). Fig. 5 shows the typical computed dependences 
of the mean mutual information (P 3 ) on the coupling 
strength e for RCM in the cases a = 2, 1.8, and 1.6, and 
GCM. At low coupling intensities, the mutual informa- 
tion is zero, indicating the absence of correlations in the 
turbulent phase. It starts then to increase and reaches 
a maximum. In GCM this maximum corresponds to the 
ordered phase. When the coupling intensity is further in- 
creased, the mutual information falls down -at the onset 
of the glassy phase for GCM- before it finally increases 
and reaches a stable high value in the synchronous state. 
The effect of a decreased connectivity translates into a 
shift of the phases to higher values of the coupling e. Fol- 
low enough v 1 synchronization is no longer possible (see 
also Fig. 3). 

Thus, though by direct counting of the number of clus- 
ters (Fig. 1) we could not see any ordering in RCM for 
e < e* , the present analysis based on the measure of the 
mutual information clearly shows that the networks have 
intrinsic dynamical organization also in this region. 
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FIG. 5. Mean mutual information ) as function of the 
coupling intensity e for a GCM (v — 1) of size N = 200 
(dashed line) and three randomly chosen networks of the same 
size, connectivity v = 0.8 and values of a as shown. The se- 
quences S\ were taken after discarding a transient of 10 3 steps. 
The GCM data is additionally averaged over 100 random ini- 
tial conditions, and the RCM over 25 independent graphs. 

Further characterization of different RCM phases is 
provided by two order parameters jq]. We examine all 
different N(N — l)/2 pairs of elements in the system and 
count how many of them are at a distance shorter than 
a fixed given precision S. The first order parameter r(S) 
is defined as 

r{6) = N(N — 1) £ 9 (* - I s * W - **(*)!) < 9 ) 

1 ' (i<3) 

where Q(x) is a step function, 9(x) = for x < and 
Q(x) — 1 for x > 0. The sum is taken over all possible 
ordered pairs (i < j). The second order parameter s(S) is 
given by the relative number of elements having at least 
one other element at a distance d < S. 

In globally coupled systems the synchronization pro- 
ceeeds until all the elements in the same cluster asymp- 
totically reach identical dynamical states. In this case, 
one can choose the highest available precision in the cal- 
culation of the two order parameters (in actual simula- 
tions it is limited by the computer precision) . As we have 
already noted, this absolute synchronization does not oc- 
cur in RCM. Instead, only clusters of elements having 
close dynamical states are formed here. Therefore the 
choice of the precision S becomes important when net- 
works are considered. 

We have computed the order parameters as functions 
of the effective coupling strength tjv for RCM and the 
respective GCM using varying precisions S. It was found 
that the results for GCM only weakly depended on the 
precision in a wide interval ICC 10 < 6 < 1CP 3 . The curves 
shown in Fig. 6(a) for GCM have been calculated us- 
ing 5 = 10~ 6 . We see that the order parameter r for 
GCM reaches, as should be expected, the value r = 1 in 
the fully synchronous (coherent) state at high coupling 
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intensities. However, a large relative number r of syn- 
chronous pairs is also found in this case at lower coupling 
intensities e in the ordered phase. Moreover, the second 
order parameter s in this phase is close to 1 indicating 
that almost all elements belong to one of the synchronous 
clusters. 

In contrast to globally coupled systems, dynamical 
clustering and synchronization in RCM is best resolved 
when an optimal small precision is employed. The plots 
shown for RCM in Fig. 6(b) have been therefore con- 
structed for the optimal precision 5 = 10~ 3 . We see that 
in the region e > e* , where fuzzy mutual synchronization 
is observed, both order parameters reach their maximal 
possible values, r = s = 1. Below the synchronization 
threshold, both order parameters rapidly decrease but 
then show a maximum. For small coupling intensities, 
the order parameters become very small (a somewhat 
larger initial level of s is explained by the fact that S is 
larger here than in Fig. 6(a) and therefore a small num- 
ber of pairs separated by the distance S is found already 
in the turbulent phase, for the random independent dis- 
tribution over the one-particle attractor). 
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FIG. 6. Order parameters r (dotted lines) and s (solid 
lines) as functions of the coupling intensity e for a GCM of size 
N — 250 (a) and for a randomly chosen network (the same as 
in Fig. 5) of size N = 250 and connectivity v = 0.8 (b). The 
employed precisions are &gcu = 10~ 6 and Srcm = 10~ 3 . 
Averaging over 100 random initial conditions is additionally 
performed for GCM. 

Thus, the behavior of the order parameters is again 
qualitatively similar in RCM and GCM. We can conclude 
that within a certain interval of coupling intensities, dy- 
namical clustering of elements occurs in these systems. 
The difference is that, in the case of RCM, the clusters 
are fuzzy and can therefore be identified only when a suf- 
ficiently low precision is used. One further difference is 
that for RCM the order parameters do not fall down so 
sharply immediately below the synchronization thresh- 
old and a significant number of elements still has close 
neighbours in this region. 



IV. DISTRIBUTIONS OVER PAIR DISTANCES 
AND DYNAMICAL CLUSTERS 

Additional information about the structure of different 
phases in RCM is provided by histograms of distributions 
over pair distances. Such histograms are constructed by 
counting at a given time moment the numbers of pairs 
with distances dij = \xi — Xj\ lying within subse- 
quent small intervals Ad. Fig. 7 shows these normalized 
histograms for one fixed randomly chosen network with a 
large number of elements (N = 1000) at several coupling 
intensities e. In the turbulent phase, a flat distribution 
corresponding to almost independent elements is found 
(Fig. 7(a)). When the coupling intensity is increased, 
some inhomogeneities start to develop in the distribu- 
tion (the mutual information and the order parameters 
also begin here to increase) . The peaks appearing later in 
the distributions indicate the onset of dynamical cluster- 
ing (Fig. 7(b)). In the situation shown in Fig. 7(c) the 
system has two fuzzy clusters. When the coupling inten- 
sity is further increased, the cluster structure is destroyed 
and the distribution characterized by a broad maximum 
at zero distance between elements is formed (Fig. 7(d)). 
As the coupling intensity e grows, this maximum gets in- 
creasingly narrow until the synchronous state is reached 
at e = e*. 




FIG. 7. Normalized histograms of distributions over pair 
distances dij for different coupling intensities (a) e = 0.1, (b) 
e = 0.25, (c) e = 0.35, (d) e = 0.45 for a randomly chosen 
network of size N = 1000 with connectivity v = 0.8. The 
histograms are obtained by counting the numbers of pairs 
with the distances falling inside subsequent intervals of width 
Ad = 0.01 at a fixed time moment t = 200 after the transient. 
The vertical and horizontal axes have the same scales in all 
these histograms and stand for the probability density and 
for the distance between states, respectively. 

The distributions over pair distances in Fig. 7 corre- 
spond to fixed time moments and therefore cannot tell us 
anything about dynamical properties of the clusters. To 
analyze the underlying dynamical behavior of the system, 
we have plotted in Fig. 8(a-d) the typical time evolutions 
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of the distance between two elements for the histograms 
displayed in Fig. 7(a-d). 

In the turbulent phase (Fig. 7(a) and 8(a)) the pair 
distance evolves in an irregular way and shows large vari- 
ations, as it can be expected for two independent logistic 
maps. At the beginning of the clustering phase (Fig. 7(b) 
and 8(b)), the elements tend to stay much closer (notice 
the change in the vertical scale) and weak aperiodic os- 
cillations are observed. 

For coupling intensities near the top value of the mu- 
tual information (Fig. 7(c) and 8(c)), the elements lock 
into periodic trajectories. Examining the trajectories of 
individual elements, we have seen in this case that all 
of them are now periodic, though different for different 
elements of the system. Thus, the system acquires rigid 
internal organization and falls into a state of frozen disor- 
der. The pair distance between two elements in a cluster 
shows in this case purely periodic variation. The clusters 
are rigid and no exchanges among them are observed. 
Figures 7(b,c) represent two elements belonging to the 
same cluster. The distance between intercluster pairs 
shows analogous behaviour (i.e. weak aperiodic oscil- 
lations or periodic dynamics) although then the typical 
separations are of order O(l). 

An interesting dynamical behavior is observed for 
higher coupling intensities, preceeding the synchroniza- 
tion transition (Fig. 7(d)). Now the elements alternate 
between short periods of partial synchronization and ex- 
cursions away from the incipient cluster (Fig. 8(d), note 
the increase in the vertical scale) . This form of behavior 
is in fact very reminiscent of on-off intermittency [ p^[ . 
Such intermittency can explain the origin of the broad 
shoulder in the histogram of Fig. 7(d): It is formed by 
the elements that temporarily find themselves during a 
large excursion from the central cluster. 




FIG. 8. Time evolution of the distance between a pair of 
elements in the same random network as in Fig. 7 for different 
coupling intensities (a) e = 0.1, (b) e = 0.25, (c) e = 0.35, and 
(d) e = 0.45. In cases (b) and (c) the elements belong to the 
same cluster. 



The previous analysis has been also carried out for a = 
1.8. The above described picture is also obtained for 
this other value of the parameter in the logistic map and 
for corresponding coupling strenghts e = 0.1, 0.15, 0.25 
and 0.29. This sequence is coherent with the different 
phases detected by the mean mutual information (P J ) 
and represented in Fig. 5. 

The formation of clusters with periodic dynamics has 
previously been observed in the ordered phase of glob- 
ally coupled logistic maps |J. As the coupling intensity 
is further increased, this ordered phase is replaced in this 
system by the glassy phase where the system has a large 
number of different attractors and its asymptotic dynam- 
ics strongly depend on the initial conditions. The glassy 
phase of GCM preceeds the final transition to the fully 
synchronous coherent state. 

An important result of our study is that the glassy 
behaviour was absent in the studied randomly coupled 
maps. When dynamical clustering was observed in this 
system, the cluster structure did not depend on the initial 
conditions and was completely determined by the archi- 
tecture of the underlying graph. Moreover, the phase 
of dynamical clustering is separated in RCM from the 
synchronous phase by the region of intermittent regimes. 

Our interpretation of this finding is that the quenched 
disorder introduced by randomly deleting some connec- 
tions -transforming the GCM into RCM- is to a certain 
extent equivalent to the introduction of a small amount 
of dynamical noise (either multiplicative or additive) in a 
globally coupled system |22| . We have checked this con- 
jecture by constructing the distributions over pair dis- 
tances for GCM and RCM in the glassy and intermittent 
phases, respectively. Moreover, we have also computed 
similar distributions for GCM where an additive or a mul- 
tiplicative noise have been included. The dynamical evo- 
lution of the noisy GCM is defined through 

N 

x\t + 1) = (1 - e)f{x\t)) + ± £ /(**(*)) + 



(10) 



where g{x l (t)) = Ti{t)x l {t) in the multiplicative case and 
g{x l (t)) — ri(t) in the additive case. We have used a 
small amplitude £ — 10~ 3 for the noise, and r^t) is a 
random number between — 1 and 1 . It is chosen anew for 
each map at each time step. 

Fig. 9(a) and (b) show in a logarithmic scale the his- 
tograms of distributions over pair distances in the glassy 
phase of the globally coupled logistic map under two dif- 
ferent choices of the initial conditions for the same cou- 
pling intensity. We see that the resulting distributions 
are very different. Note that both distributions have been 
averaged over time, so that the peaks and irregularities in 
these figures reveal the persistent structure of the under- 
lying attractors. Fig. 9(c) shows how these distributions 
are influenced by introducing into the GCM a weak ad- 
ditive (dashed line) or multiplicative (bold line) noise of 



7 



intensity ( = 1CP 3 according to Eq. (7). The noises 
wash out the fine jagged structure of the distribution 
and, more importantly, make it independent of the ini- 
tial conditions. The resulting distributions become then 
clearly similar to the respective distribution we obtain at 
the same coupling intensity for RCM, Fig. 9(d). 




0.0 0.2 0.4 0.6 0.8 1.0 

d 



FIG. 9. Normalized histograms over pair distances for a 
GCM of size N = 50 and two different initial conditions (a,b), 
for the same GCM in the presence of an additive (dashed line) 
or multiplicative (bold line) noise of intensity £ = 1CP 3 (c), 
and for a randomly chosen network of the same size with 
connectivity v = 0.8 (d). The coupling intensity is e = 0.45 
in all these plots; time averaging is additionally performed. 

We would like to emphasize that this parallelism be- 
tween quenched disorder and dynamical noise might hold 
in the intermittent phase, but of course the mechanisms 
leading to the dynamical behaviour observed in other 
phases cannot be (at least solely) ascribed to noise. In the 
clustering phase, for instance, it should be clear that the 
fixed (though disordered) structure of the network plays 
the main role. In this sense, we have analysed the dy- 
namical behaviour of the clusters by measuring in a fixed 
network how many elements belong to a certain cluster 
and how many clusters are formed at each time step. A 
cluster C m (S;t) for a given precision 6 and at time t is 
formed by k m elements, m — 1, . . . , M such that all of 
them have at least another element of the cluster at a 
distance d < S, that is d = \xi(t) — Xj(t)\ < S, Vi and 
some j G C m in order to say that also i G C m . 

In Fig. 10 we represent the size of all clusters in a 
network with N — 50 elements as a function of time (for 
a fixed precision 5 = 0.1). Four different values of e in 
the clustering phase have been chosen. In the first plot 
(Fig. 10(a), e = 0.25), the elements tend to cluster but 
the groups are still relatively unstable. A closer inspec- 
tion of the clustering dynamics reveals that a group of 
19 elements keeps stable in time, while another cluster 
containing the 31 remaining elements splits often in sub- 
groups of sizes (17,14), (18,13) or (20,11) among others. 
If the coupling strength is increased, also the stability of 



the clusters increases, and their average lifetime becomes 
longer. In Fig. 10(b), for e = 0.28 we observe that, in 
fact, larger and more stable clusters are formed. Now 
the elements are divided into a cluster with 20 elements 
and a second group with 30 that often splits into two 
subclusters with 27 and 3 elements, respectively. Some 
irregularities in the dynamics are also found. In the hard- 
locking phase for this network (Fig. 10(c), e = 0.30) two 
stable clusters with 24 and 26 elements are formed. As 
already discussed, the maps display periodic trajectories 
in this narrow parameter region. For a slightly larger 
e = 0.32, in Fig. 10(d), we see a first stable group with 
27 elements and a second one with 23 including a weakly 
coupled map (which periodically leaves the cluster). 



(a) - 




(b) 
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FIG. 10. Size and stability of synchronous groups in the 
clustering phase. The dynamical behaviour of a fixed network 
with N = 50 elements and v = 0.8 is analysed for coupling 
strengths e = 0.25 (a), e = 0.28 (b), e = 0.3 (c), and e = 0.32 
(d). The precision is in this case 8 — 0.1. Weakly coupled 
elements coexist with long-lived clusters which often split into 
smaller subgroups. 

We have observed that the transient time required for 
the system to fall into the clustering phase increases with 
the increase of e. At the same time, it is seen that ele- 
ments have a trend to condense in a single group, as 
revealed by the presence of some time steps where the 
cluster size equals N. This is not found to happen at the 
beginning of the clustering phase, for values of e < 0.27. 
Finally, for e 0.35 (in the network of Fig. 10) the maps 
do not form clusters any more and the intermittent phase 
begins. 

V. PARTITIONS INTO DYNAMICAL CLUSTERS 

In this section we more closely examine the struc- 
ture of the dynamical clustering phase in RCM. In this 
phase, groups of maps moving together in a robust way 
and forming long-lived clusters have been observed. The 
emerging cluster structure is biased by the connection 
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patterns of the underlying network. To demonstrate this, 
we introduce relative connectivities that are defined be- 
low. 

Let us suppose that at time t and with precision S 
our system separates into M clusters, C m (5;t), m 

1 M, each of which containing k m elements. The 

relative connectivity inside a cluster m is then defined as 



1 



'«// 



ukm (fori 



i 



i,jeC m , i^j (11) 



where the sum is taken over all pairs of elements belong- 
ing to this cluster. The relative connectivity between two 
different clusters / and n is given by 



Ln 



1 



e ff vhk. 



■ 2^ , 



i e Ci, j e C n 



(12) 



Thus defined, the relative connectivities are equal to 1 if 
the characteristic connectivity inside a cluster or between 
two clusters are exactly the same as the average con- 
nectivity v of the entire network. Positive deviations of 
the connectivity inside a cluster (z/J . > 1) indicate that 
this cluster contains elements which are more strongly 
connected than on the average. Respectively, when the 
relative connectivity between two clusters is decreased 
(Veff < 1)7 this shows that the elements belonging to 
these two separate clusters are less connected than on 
the average and vice versa. 

We can also define the average intercluster relative con- 
nectivity of the entire network 



Inter 
y eff 



1 



M(M - 1 



M 

El.n 
V eff 

n,l—l;n^l 



(13) 



and its average intracluster relative connectivity 



, Intra 
V eff 



M 



M 



(14) 



To characterize the cluster structure of the partially 
ordered phase, we fix the coupling intensity e and con- 
sider the state of the whole system at a given time mo- 
ment. By varying the precision <5, we obtain a hierarchy 
of cluster partitions, as seen with different resolutions. 
For each resolution level, its relative connectivities are 
then calculated. Fig. 11 presents the emerging hierarchi- 
cal structure of dynamical cluster partitions for a system 
of N = 50 elements with coupling strength e = 0.23 at 
three different precisions S. The numbers between brack- 
ets correspond to the number of elements in each partic- 
ular cluster. The numbers inside the clusters are their 
relative connectivities and the numbers on the links be- 
tween the clusters yield the relative connectivity between 
them. Though this pattern refers to a particular time 
moment (after a long transient), it remains fairly stable 
in time. We see that, as the precision S is improved, the 



clusters split into smaller subclusters, revealing a hierar- 
chical tree-like structure It can also be observed in 
Fig. 11 that the relative connectivities inside a cluster 
exceed one, whereas the relative connectivities between 
the clusters are typically smaller than one. This indi- 
cates that the partition into dynamical clusters is biased 
by the pattern of connections in the underlying network, 
i.e. the elements belonging to a same cluster would gen- 
erally have more connections inside this cluster than with 
the elements belonging to other dynamical clusters. 



5=0.01 




FIG. 11. Hierarchical structure of dynamical clusters for a 
randomly chosen network of size N = 50 with the connectivity 
v = 0.8 at e = 0.23. 

To check more accurately this suggestion, we have cal- 
culated inter- and intracluster connectivities v I e ^ ra and 
i/inter £ or a i ar g er system with N = 250 and v = 0.8 
at e = 0.3 with the precision 6 = 0.1. These properties 
were averaged over 10 4 independent graphs. The aver- 
age intracluster connectivity was in that case ^f// ra ^) 
1 .01 .')! 1 : . that is, slightly higher than numerically gener- 
ated average connectivity, (v e ff) = 1-0000(1). The av- 
erage intercluster connectivity was ^e// er ^) = 0.987(1) 

and thus lying below (v e ff)- Figure 12 shows the nor- 
malized probability distribution over intracluster (solid 
line) and intercluster (dashed line) connectivities in the 
studied ensemble of 10 4 independently generated graphs. 
The maxima of the two distributions are slightly shifted. 
But, perhaps even more important, we see that the dis- 
tribution of the intercluster connectivities is significantly 







broader and has a wide shoulder extending towards lower 
connectivities. 
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FIG. 12. Statistical distributions over intracluster (solid 
line) and intercluster (dashed line) effective connectivities in 
the clustering partitions in an ensemble of 10 4 independently 
generated random networks of size N = 250 with mean con- 
nectivity v = 0.8 and coupling strength e = 0.3. In this case, 
we have taken the precision 5 = 0.1. 

VI. DISCUSSION 

Our numerical analysis reveals that mutual synchro- 
nization and dynamical clustering represent a typical and 
robust form of collective dynamics in random networks 
of coupled chaotic elements. The mutual synchronization 
remains possible when almost half of all potential connec- 
tions between the elements are deleted and the dynam- 
ical clustering may be found even for the more sparsely 
connected networks. For very low connectivities, we have 
seen that the system could no longer synchronize and the 
dynamical clusters in the partially condensed phase be- 
came less stable, i.e. their lifetimes were getting shorter. 

The different dynamical phases of RCM have been 
specified and compared with their counterparts in GCM. 
The essential differences have been noticed in the cluster 
structure and in the dynamics, as well as in the depen- 
dence on the initial conditions. A rich clustering struc- 
ture, depending on the network architecture, was ob- 
served in RCM. The analog of the glassy phase of GCM 
was however not found in the investigated randomly cou- 
pled maps, i.e. we have not seen that the final attractor 
depended on the initial conditions for any set of parame- 
ters. In this sense, the quenched disorder of the random 
network appears to play a role similar to that of noise 
in this phase. It might also be that the transition from 
synchronization to the intermittent phase in RCM would 
admit a characterization in terms of a blowout bifurca- 
tion |l4| , and that the intermittent phase that we observe 
immediately after the synchronous state be in fact a case 
of on-off intcrmittcncy. This picture would be consistent 
with our numerical results and with the fact that the 



disorder in the network architecture destroys the degen- 
eracy of the dynamical matrix in GCM || and generates 
a whole hierarchy of Lyapunov exponents. 

Though our investigations were made only for networks 
formed by coupled logistic maps, similar results would 
probably hold for networks made of other chaotic maps 
or elements with continuous chaotic dynamics. Indeed, 
the behavior in globally coupled logistic maps strongly 
resembles what is found in various globally coupled pop- 
ulations of chaotic dynamical systems [^|j8| . 

The systematic study of RCM implies the analysis of 
the behavior of the system under the change of four rel- 
evant parameters: The average connectivity v, the cou- 
pling strength e, the parameter of the individual map 
a, and the system size N . Our main interest in this 
study was to introduce Randomly Coupled Maps and to 
give some insights on the role of the network architecture 
in the dynamics. Hence, we have mainly investigated 
the two parameters v and N, and reanalysed the known 
phases for GCM when e varies from zero to unity. Many 
of our investigations were performed with control param- 
eter a = 2 of the logistic map. This value is somewhat 
special, since at a > 2 the trajectories become infinite 
and the chaotic attractor disappears in a boundary crisis 
. Other simulations for smaller values of the control 
parameter a show a similar qualitative behaviour. 

In our study, the networks were generated by inde- 
pendently choosing with a certain fixed probability the 
connection between any two elements. Thus constructed, 
the connection patterns are random, but statistically uni- 
form. We have analysed systems with sizes ranging from 
a few elements to N — 2 11 . Generally, the synchroniza- 
tion threshold depended not only on the system size and 
on the average connectivity, but also on the particular 
architecture of a chosen network. We have seen, how- 
ever, that variations in the synchronization threshold for 
networks with the same mean connectivity became much 
weaker when the network size was increased, and that 
the distance to the mean field threshold given by GCM 
had the functional form (e* — Cqcm) — N^ 1 / 2 . The ex- 
istence of a universal synchronization threshold in such 
randomly generated networks in the limit vN — > oo is 
thus expected. The statistical uniformity, introduced in 
this paper through the independent choice of individual 
connections, is a special feature that should not neces- 
sarily be present in complex networks. Natural networks 
may have various topological structure s [P5| which can 
also result from evolutionary processes |26|J27|] . It would 
be interesting to see how synchronization and the dynam- 
ical clustering phenomena arc influenced by such struc- 
tures. 

We have found that the network architecture biases 
the partitions of the network into dynamical clusters and 
determines interactions between the clusters which lead 
to their collective dynamics. This puts forward the task 
of engineering the networks with the desired dynamical 
clustering properties. One can apparently design sys- 
tems that would display an arbitrarily chosen partition 
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into several exactly synchronized clusters (see the Ap- 
pendix) . The collective dynamical behavior can represent 
an important practical function of a network. The evo- 
lution of a network, proceeding through random muta- 
tions, may then be guided towards the optimization of its 
collective dynamics. Indeed, examples of dynamical net- 
works that evolve to reproduce given temporal 'melodies' 
have already been constructed Pq^9(| . We want there- 
fore to emphasize that the evolution of networks can also 
be steered to reach better synchronization properties or 
to approach a certain dynamical clustering structure. 

Finally, we note that when the dynamical clustering 
is taking place, coherent clusters can be interpreted as 
super elements that form an emerging dynamical net- 
work of a higher structural level. Taking into account 
the large variety of clustering partitions and their sensi- 
tivity to the coupling intensity, RCM systems may thus 
be viewed as a living space that supports different dy- 
namical (meta)networks and may retrieve a particular 
such network under appropriate changes of the control 
parameters. 
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APPENDIX 

It was earlier noted that the network architecture in- 
fluences the critical coupling intensity e* at which syn- 
chronization first appears and favours certain preferred 
partitions of elements into dynamical clusters. In this 
Appendix we analyze the role of the network architec- 
ture in the dynamical clustering phenomena for small 
networks consisting of only N — 5 elements. If the net- 
work connectivity is fixed at v = 0.6, there are just four 
such networks shown in Fig. 13. We have systematically 
investigated their synchronization properties for various 
values of the control parameter a of the logistic map in 
the interval from 1.42 to 2 with increment Aa = 0.02 
and for the coupling intensity e in the interval from 
to v with increment Ae = 0.01. In this case we have 
considered that two maps are synchronized if they have 
exactly the same state. This is now licit because of the 
high degree of symmetry of the networks. The percent of 
parameter pairs leading to each of the possible clustering 
configurations is displayed in Table 1. In Fig. 13, ele- 
ments with the same symbol synchronize (i.e. they form 
a cluster) with the higher probability. Different symbols 
stand for different clusters. The value of e at which the 
elements in each of the networks synchronize can take 
a wide spectrum of values. For instance, for a = 1.6 it 



changes from 0.56 (case(A)) to 0.96 (case(B)) [p0[ . More- 
over, we have found that e* is a non-monotonic function 
of a, and can be even decreasing depending on the graph. 
Hence, if N is small, each network has to be treated in- 
dependently (as in the example here analyzed). 

We see that indeed in the majority of cases the clus- 
tering partition follows the pattern of connections in the 
graph. Moreover, some of the potential highly asymmet- 
rical partitions have never been observed (like the par- 
tition into (1,5) (2,3) (4) for the graph A). This shows 
that the symmetry of connections inside a graph plays an 
important role in the dynamical clustering phenomena. 

Looking at the graph A, we see that the dynamical 
equations of its elements are not changed under the re- 
labeling (1,2,3,4,5) — * (2,1,4,3,5), reflecting the sym- 
metry with respect to permutations in the matrix Ty for 
this graph. It seems highly plausible that synchronous 
clusters would generally be much easily formed by indis- 
tinguishable elements, defined as those whose dynamical 
equations are identical under permutations. This is the 
reason that allows synchronization to be of the hard lock- 
ing (exact) type in this case. The numerical results shown 
in Table I support this statement. 

Consider, for example, the situation in which the clus- 
ters (1,2), (3,4) and (5) have been formed in (A). Let us 
call x x {t) = x 2 (t) ee x, x 3 (t) = x 4 (t) = y and x 5 (t) ee z 
and write down the dynamic equations for the new 3- 
cluster system: 

x(t + l)= (l-^) /W + ^e/(y) 



.'/('+ I I') f(y) + ^e(f(x)+f(z)) 



z(t + l) 



6 



f(z) + -ef(y) (Al) 



On the one hand, we do not observe any further (at 
least trivial) symmetry in these equations. On the other 
hand, the problem of global synchronization in our orig- 
inal graph has moved to the problem of synchronization 
of asymmetrically connected oscillators which moreover 
have different values for the coupling strength. Accord- 
ing to this result, we believe that the problem of syn- 
chronization and clustering in RCM might also include 
asymmetrically connected networks and probably some 
cases of coupled systems with a distribution n(e) of e 
values p|jl|. 
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FIG. 13. Four different possible configurations of a network 
with N = 5 elements and z/ = 0.6. In each graph, favoured 
synchronizations are displayed using the same symbol. 
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TABLE I. Percent of realizations (%) for the clusterings shown and for each of the networks in Fig. 13. 



Network 


Clusters 


% 


Network 


Clusters 


% 


A 


(1,2,3,4,5) 


5 


B 


(1,2,3,4,5) 


1.5 




(1,2,5) (3,4) 


2 




(1,2,4) (3) (5) 


1 




(1,2) (3,4,5) 


1 




(1,2) (3) (4) (5) 


75 




(1,2) (3,4) (5) 


59 




(1,2,3,4) (5) 


1 




turbulent 


33 




turbulent 


21.5 


C 


(1,2,3,4,5) 


2 


D 


(1,2,3,4,5) 


4 




(1,2,3,4) (5) 


1 




(1,2,3,4) (5) 


1 




(1,3,4) (2) (5) 


17 




(1,2) (3,4) (5) 


57 




(1,3,4) (2,5) 


0.5 




(1) (2) (3,4) (5) 


1.5 




(1,2) (3,4) (5) 


0.5 




(1,2) (3) (4) (5) 


1.5 




(1) (2) (3,4) (5) 


17 




turbulent 


35 




turbulent 


62 
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